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Abstract 


The aim of this master thesis is to investigate the influence of intergalactic dust 
on astronomical observations. The current cosmological paradigm of an accelerating 
cosmic expansion is supported by observations of Type Ia supernovae. However, the 
light emitted by these and other cosmological sources is not only redshifted by cosmic 
expansion but will also interact with matter along the light path which affects the light. 
Especially intergalactic dust can lead to additional reddening and dimming of distant 
sources due to light scattering or absorption. This yields systematic contaminations to 
cosmological measurements. Using quasars as backlight sources we seek to estimate the 
attenuation of their observed spectra caused by intergalactic dust in order to determine 
the existence and properties of the intergalactic dust. In this thesis we will analyse 
the methods and tools used as well as discussing the dust properties and how they 
affect the results. This thesis shows that intergalactic dust exist, and has a notable 
effect on cosmological measurements and favors dust grains consisting of silicate with 
a size around 0.1 wm with Qaust around 5 — 7 x 10~°. The conclusion is in agreement 
with previous studies of the same matter. With the increasing precision of cosmological 
measurements systematic effects such as intergalactic dust is a factor that has to be 
taken into account. 


Sammanfattning 


Syftet med denna masteruppsats ar att undersdka de effekter intergalaktiskt stoft 
har pa astronomiska observationer. Den ledande kosmologiska moddellen som sager att 
universum expanderar ar st6ttad av observationer av Typ la supernovor. Ljuset som 
kommer fran supernovorna ar déremot inte bara paverkade av den kosmiska expansio- 
nen men kommer ocksa interagera med materia pa vagen som paverkar ljuset. Speciellt 
intergalaktiskt stoft kan leda till ytterligare dimning av ljuset och att ljuset blir rddare. 
Detta ar pa grund av spridningen eller absorptionen av ljuset nar ljuset interagerar med 
stoftpartiklarna. Detta leder till systematiska fel i kosmologiska matningar. Genom 
att anvanda kvasarer som en ljuskalla forséker vi uppskatta ljusd&émpningen pa deras 
spektrum som orsakats av intergalaktiskt stoft. Detta gors for att bestamma ifall in- 
tergalaktiskt stoft existerar och vad dess egenskaper dr. I denna uppsats gjordes en 
analys av metoden, beraékningsverktygen som anvandes och nagra egenskaper hos stoft. 
Slutsatsen av denna analys visar pa att intergalaktiskt stoft existerar och har en mat- 
bar effekt pa kosmologiska observationer. Resultatet visar att modellen f6redrar stoft 
av silikat med storlekar runt 0.1 wm med Qaus: runt 5 — 7 x 10~°© vilket st&émmer dv- 
erens med tidigare studier. Med den dkade pressionen av kosmologiska matningar blir 
systematiska effekter, som den fran intergalaktiskt stoft, Annu viktigare att ta i hansyn. 
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1 Introduction 


Cosmological observations suggests that the universe is expanding at an accelerated rate 
(Hubble 1929). This expansion is supported by observations of Type la supernovae which 
have well-defined luminosity properties and lightcurves resulting in that they can be used 
as standard candles (Riess et al. 2016). In an expanding universe the lightcurves produced 
by these standard candles would be shifted towards the red part of the spectrum due to the 
light travelling through expanding space, resulting in a cosmological redshift. The further 
away they are the further shifted towards the red the lightcurves are since the space the light 
travels trough has underwent greater expansion. However there are other factors which could 
induce a redshift of the lightcurves. One of these factors are the matter the light interact 
with along the path to us, especially intergalactic dust can lead to additional reddening and 
dimming of distant sources due to the light being scattered or absorbed and re-emitted by 
the dust grains (Johansson 2015). The interaction between the light and the dust grains thus 
yield a systematic contamination in cosmological measurements. In Figure 1 we can see how 
cosmic dust reddens the light from sources behind the comic dust. This contamination could 
contribute to the discrepancy between differently attained cosmological constants (Mortsell 
et al. 2021). 

Cosmic dust can spawn from various sources. These sources include Type IT supernovae, 
Asymptotic Giant Branch (AGB) stars, galaxy and star formation, etc. However for the 
dust to dilute the intergalactic media the dust needs to leave its origins and the galaxies 
from which they were born. To be able to do this the dust must have acquired sufficient 
kinetic energies in their making such that they can escape the potential well of the galaxy. 
This means that supernovae are potential producers of intergalactic dust since the dust 
created in such events can be ejected with great kinetic energies. With time more and more 
supernovae and other astronomical processes which produces dust will occur thus enriching 
the intergalactic medium (IGM) with dust over time (Johansson and Mértsell 2012). 

Cosmic dust is also a key ingredient in understanding the formation and evolution of 
stars, since stars are formed in regions full with dust as well as cosmic dust being important 
for many evolutionary steps of stars (Karttunen H. et al. 2017). Dust also plays a role in 
the formation of galaxies and their evolution. So in order to gain better understanding of 
both stars and galaxies knowledge about cosmic dust plays a big part. Lets summarize why 
cosmic dust is important: 


e It plays a big role in the formation and evolution of stars. 
e It is key to understanding galaxy formation and evolution. 
e It causes a reddening of sources behind the dust. 


e As the precision of cosmological measurements increases systematic effects, such as 
those from intergalactic dust, plays a larger role. 


e It could help to understand the gap in the Hubble tension. 


In this thesis we investigate the influence of intergalactic dust on observations, as well 
as constraining the properties of intergalactic dust such as the dust grain size. We do this 
by constructing a forward model where we model a simulated quasar spectrum after which 
we apply our intergalactic dust attenuation model onto the spectrum. Once we can generate 
intergalactic dust attenuated quasar spectra we compare these to SDSS data of observed 
quasars by comparing the colours of the simulations and the observed data. We then fit the 
input parameters of the dust model to the SDSS data by doing a statistical analysis of several 
different dust models with varying parameters. This thesis assumed further simplifications 





Figure 1: A Bok globule which is a dense region of cosmic dust and gas where star formation 
may take place. As seen in the image the cosmic dust cloud reddens the light from the 
sources behind the dust cloud. Image source: Bok globule B68, ESO 2001 


of the model that was initially planned not to be included, more on this in Section 8. Parts 
of this thesis is an updated version of previous work (Kjellqvist 2021). 

Let me end the introduction with the inspirational quote from Carl Sagan: 

“Even through your hardest days, remember we are all made of stardust.” - Carl Sagan 


2 Stellar evolution 


When investigating intergalactic dust and its influence on observations, especially over dif- 
ferent redshift, it is necessary to take the production of dust into consideration and how that 
has changed the number density of dust grains over time. The creation of intergalactic dust 
grains are mainly linked to different evolutionary stages of stars (Masataka Fukugita 2011, 
Gall, Hjorth, and Andersen 2011), so lets take a look at stellar evolution. 

All stars are formed from collapsing clouds of gas and dust. These massive clouds are often 
called molecular clouds or nebulae. These protostars contract over the course of millions of 
years until they finally enters a state of equilibrium and enters the so called Main-Sequence 
(MS). How long a star spends on the MS depends on the mass of the star, where more 
massive stars are shorter lived compared to less massive. In Table 1 lifespans for an array of 
masses (Bertulani 2013). 





Figure 2: The evolutionary tracks of stars along the MS shown in a Hertsprung-Russel (HR) 
diagram. The evolution of the stars depends on where they start on the MS. In the plot sun 
is seen on the MS in the yellow part. The giant and supergiant stars are also seen in the 
diagram, as well as the white dwarfs. Image source: Main Sequence, Universe Today 2010. 















































Mass (solar masses Mo) Time [years| Spectral type 
60 3 million O03 
30 11 million O7 
10 32 million B4 
3 370 million A5 
1.5 3 billion F5 
1 10 billion G2 (Sun) 
0.1 1000s billions M7 
Table 1: The lifetime of stars with different masses along with their spectral type (Bertulani 





2013). The spectral type is a classification based on spectral characteristics of the star, which 
depends on the temperature of the star, O being the coolest and M the hottest. 


Nuclear fusion is what powers the stars during the majority of their existence. Initially 
this is via fusion of hydrogen into helium in the core of MS stars. As the hydrogen is 
consumed in the core the star will start to move away from the MS. How the star moves 
through and beyond the MS depends on the mass of the star. One can thus divide the stars 
into low (<0.6 Me (solar mass)), mid (0.6-8 Me) and massive (>8 Mo) sized (Iben 1967). 
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Figure 3: The entire stellar life cycle for different mass stars, starting from the birth in the 
molecular clouds and ending with the remnants. Image source: Star Life Cycle Chart, R.N. 
Bailey 2017 


2.1 Low-mass stars 


Due to the extremely long lifespan of low-mass stars we have not observed what happens to 
these stars once they stop producing energy through fusion. Since the universe is around 13.8 
billion years old and the lifespan of these stars are sometimes several orders of magnitudes 
larger we can only investigate the evolution of these stars through theoretical models and 
computer simulations. Models of red dwarfs with a mass of 0.1 Mo suggest that the such 
stars will stay on the MS for around six to twelve trillion years. During this period they 
would gradually increase in both temperature and luminosity, and finally moving of the MS 
taking several hundred billion years to collapse into a white dwarf, thus never evolving into 
ared giant (F. Adams, P. Bodenheimer, and G. Laughlin 2005). 

Stars which are slightly more massive will eventually expand into red giants but their 
cores will not become massive enough to reach the temperatures required for burning of 
helium thus never reaching the tip of the red giant branch (see Figure 4). These stars will 
also eventually collapse into white dwarfs (Gregory Laughlin, Peter Bodenheimer, and F. C. 
Adams 1997). 


2.2. Miid-sized stars 


Mid-sized stars go through several stages during their evolutionary track. Figure 4 shows 
their evolutionary track through the HR diagram along with a brief explanation of each 
phase. 
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Evolution of a 1 Mo star 
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Figure 4: The evolutionary track of a 1 Mo star shown in a HR diagram. It starts at the main 
sequence moving into the subgiant branch, then into the red giant branch eventually reaching 
the Asymptotic giant branch (AGB), and lastly reaching the post ABG phase and ending as 
a planetary nebula and a white dwarf. Image source: Evolutionary track, Lithopsian 2016 


2.2.1 Subgiant phase 


Once the hydrogen in the core of MS stars are depleted the stars leave the MS and enters 
the subgiant phase. In this phase the star begins to fuse hydrogen in a shell around its 
helium core, continuously increasing the mass of the core as more helium is created in the 
shell. During this period the star expands and cools. This process occurs over a period of 
several million to one or two billion years depending on the mass of the helium core. As the 
convective envelope of the star expands it will eventually begin to increase in temperature 
and luminosity and expand into the red giant branch (Ryan and Norton 2010). 

At acertain point the helium core becomes sufficiently massive such that fusion of helium 
begins. 


2.2.2 Red giant & the horizontal branch 


How the star evolves in the red giant phase is highly dependent on the mass of the star. For 
stars with masses < 2.3 Mo their cores will continue to grow and increase in density to a 
point where it becomes degenerate, meaning that the electron degeneracy pressure becomes 
significant due to the Pauli exclusion principle. This will lead to a uniform temperature 
increase, due to the high conductivity of the degenerate gas, eventually reaching 100 million 
degrees which is sufficient for helium burning to start via the triple alpha process (Equations 
1 and 2). The sudden start of helium burning will occur simultaneously throughout the entire 
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core which will lead to a great increase in temperature. Due to the core being degenerate 
it will not expand leading to a greater increase in temperature which only accelerates the 
fusion of helium. Eventually the temperature reaches a point where the electron degeneracy 
pressure is not required to keep the core from collapsing, which removes the degeneracy of 
the gas. This causes the core to expand rapidly and only a couple of seconds after the start of 
helium burning there is an explosion, the so called helium flash (Karttunen H. et al. 2017). 


$He +4 He § Be (1) 
®Be+5 He oP C+ 2y (2) 


This might sound flashy (is puns allowed in a master thesis’), but the energy released in 
the explosion is absorbed by the outer layers and turned into potential energy of the core, 
which eventually leads to a non-degenerate core which burns helium into carbon. After the 
helium flash the star reaches the horizontal branch, which gets its name from the fact that 
the luminosity does not change much while the temperature may vary (Karttunen H. et al. 
2017). 

For stars with a mass in the 2.3-8 Mo range the temperature of the cores is greater than 
the temperature of the smaller star’s cores. They also have less dense cores, which results in 
the core never entering a stage of degeneracy. This means the cores can start their helium 
burning in a non-catastrophic way which does not result in a helium flash (Karttunen H. et 
al. 2017). 


2.2.3 Asymptotic giant branch (AGB) phase 


Once the helium is depleted in the core the core will consist of carbon and oxygen. At 
this point a helium layer will form outside of this core where helium burning will continue. 
Outside of this helium layer a layer of hydrogen exists where hydrogen burning continues. 
During this stage the star will move in a similar track as it did during the red giant phase, 
but at higher temperatures, thus it is called the asymptotic giant branch. During this phase 
the energy production of the star is dominated by either the hydrogen or the helium shell 
in an alternating fashion. Once sufficient helium has been created in the hydrogen shell it 
will drop down to the helium shell igniting more rapid burning of helium. A star with two 
different layers of fusion is an unstable one and the star enters a thermally pulsing phase 
where material will periodically be ejected into space. This leads to the creation of a lot of 
dust (Karttunen H. et al. 2017, Sackmann, Boothroyd, and Kraemer 1993, Gall, Hjorth, and 
Andersen 2011). 


2.2.4 Post-AGB phase 


These pulsations and stellar winds will continue until all of the stars outer layers are expelled 
as a planetary nebula. These planetary nebulae are rich of heavy elements such as carbon 
and oxygen, and as it moves away from the star it cools allowing the creation of dust particles 
and molecules to take place. The remaining star never become hot enough to star carbon 
and oxygen burning, resulting in the star cooling down and becomes a white dwarf consisting 
of carbon and oxygen (Karttunen H. et al. 2017, Gall, Hjorth, and Andersen 2011). 


2.3. Massive stars 


There are some similarities between the evolution of intermediate massed stars and massive 
stars, but there are some very clear differences. Generally the more massive a star is the 
faster all of the evolutionary processes occur and the shorter it is lived. This shorter life- 
span leads to massive stars being more difficult to find, and is usually found in young galactic 


12 


structures. Massive stars starts, as for any other star, with hydrogen burning in their core. 
As that is depleted the star leaves the MS and starts to expand into a supergiant. They also 
start to burn helium in their core in a non-catastrophic manner in contrast to lower massed 
stars. They also differ from lower massed stars in that they are able to fuse elements beyond 
helium, meaning that they do not eject their outer layers through shell burning of hydrogen 
and helium. Instead they continuously fuse heavier and heavier elements starting with carbon 
and then oxygen and silicon burning, creating layers after each burning process is completed. 
In each layer the burning of that element continues. The most massive stars, larger than 15 
Mo, are able to fuse elements all the way to °°Fe creating onion-like structures with several 
shells with different material where shell burning occurs. These are not stable structures and 
will eventually lead to a collapse of the star when the core exceeds the Chandrasekhar limit 
of about 1.4 Mo. At this point the electron degeneracy pressure is not enough to overcome 
the gravitational pull (Karttunen H. et al. 2017). 

During the collapse the core will collapse on a faster time scale compared to the outer 
layers. The collapse results in a great increase in temperature in the outer layers where 
unburnt nuclear fuel resides. This leads to a runaway fusion reaction resulting in explosive 
burning releasing great amounts of energy in a core collapse supernova supernovae. The 
majority of the energy are in the form of neutrinos, but great amounts of matter is also 
ejected leading to the creation of dust (Karttunen H. et al. 2017, Gall, Hjorth, and Andersen 
2011). 

The very dense collapsing core will finally either turn in to a neutron star, a star consisting 
almost entirely of neutrons in a degenerate state due to the extreme pressure, or a black hole 
depending on if the mass exceeds the Tolman-Oppenheimer-Volkoff limit (Karttunen H. et 
al. 2017). 


3 Our expanding universe 


Since we will investigate the influence of intergalactic dust on a cosmological scale we have 
to take the expansion of the universe into account. In 1929 Hubble detected that far away 
galaxies seemed to be moving away from us with a recession velocity (Hubble 1929). This 
was noted due to the galaxies appearing redder than they should, implying that a Doppler 
shift towards the red part of the spectrum has occurred. A solution to these velocities were 
proposed saying that the universe is expanding. This resulted in Hubble’s law, v = Hod, 
where v is the recession velocity, Hp is Hubble’s constant and d is the distance to the light 
source. Hubble’s constant are in units kms~! Mpc~! and measures the expansion rate of the 
universe as a function of distance (Riess et al. 2016). 

For nearby objects their redshift may originate from peculiar velocities, like velocities 
due to gravitational influences. In contrast to redshift, where objects move away from us, 
sometimes the objects may be moving towards us leading to a blueshift (Barbara Ryden 
2006a). However at larger distances the redshift occurs due to the expansion of space in 
between us, the observer, and the object we are observing. This cosmological redshift depends 
on how much the space has expanded since the light we observed was emitted, and can be 
calculated with 

a(to) 


- — 
ar (3) 





where z is the cosmological redshift, a(t) is the cosmological scale factor which describes 
how much the universe has expanded as a function of time. t is the time when the light was 
emitted and to is the time when it was observed. One can use the cosmological scaling factor 
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to describe the Hubble parameter H(t) (Barbara Ryden 2006a), 


_ a(t) 

H(t) = (4) 
where a(t) is the time derivative of the scaling factor. Calculations of data from Cepheid 
variables and other astrophysical sources have determined the local value of the Hubble 
constant to be 73.24 + 1.74kms~'Mpc~! (Johansson 2015, Riess et al. 2016). There are 
however a second method of calculating the Hubble constant using the Planck satellite and 
inferring the constant from lensing measurements of the Cosmic Microwave Background. The 
interesting part is that this second method infers a value of about Hyp = 67kms~!Mpc7!, 
hence differing from the inferred value from Cepheids. The tension between the two different 
measurements of the same "constant" is a known problem, or mystery, within cosmology 
since it still is unclear why the measurements differ (Wu et al. 2020). It could be a sign 
of new physics but it could also be a sign of so far unknown systematic effects in the data 
analysis. Recent studies suggest that inaccurate colour-luminosity calibration caused by 
extinction from dust denies the precision required to claim a substantial tension between the 
two methods Mortsell et al. 2021. 





3.1 The cosmological model 


Our cosmological model builds upon Einstein’s theory of general relativity which describes 
the 4-dimensional space time structure of the universe and how it evolves. This evolution 
depends on the matter distribution of the universe, and the interaction between time and 
matter (i.e space-time) is explained in Einstein’s field equations (Einstein 1916). One can 
apply various cosmological models for different universes to the field equations and it would 
result in equations that would describe how that universe evolves as well as how matter and 
time would interact with each other. The cosmological model that best fits our universe 
postulates the so called cosmological principle. The cosmological principle states that the 
spatial distribution of matter in the universe, on large scales, is both isotropic as well as 
homogeneous. This means that on large scales the universe looks the same at every location 
and at every angle. Astronomer William Keel stated: 


"The cosmological principle is usually stated formally as "Viewed on a sufficiently 
large scale, the properties of the universe are the same for all observers.’ This 
amounts to the strongly philosophical statement that the part of the universe 
which we can see is a fair sample, and that the same physical laws apply through- 
out. In essence, this in a sense says that the universe is knowable and is playing 
fair with scientists." - William C. Keel 2007 


The Friedmann-Lemaitre-Robertson-Walker (FLRW) model assumes the cosmological 
principle for a simple matter distribution. Applying the FLRW model to Einstein’s field 
equations results in the Friedmann equations which describes the dynamics of the universe 
(Friedmann 1922), 








= (2) =~ (o+ 30 (6) 


where H is the Hubble parameter, a the cosmological scale factor, G is the gravitational 
constant, p and p is the total energy density and total pressure of the Universe respectively, 
and k is the curvature of the Universe. 
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Here it often is useful to divide the total energy density of the Universe into its con- 
tributing parts, p = Pm +pa+t pr. Here m signifies matter, which can be divided into dark 
and baryonic matter, A signify dark energy, and r which signify relativistic matter, meaning 
radiation. An equation of state parameter can be introduced for each of these parameters, 
w = p/p where 


p= poa *), (7) 
For matter (both baryonic and dark matter) the energy density is zero resulting in w,, = 0, 
for radiation w, = 1/3 and for dark energy w, = —1. Normalizing this density to the critical 


density of the universe we can now introduce the density parameter Q = p/Perit, where 
Perit = 3H?/87nG is the critical density of the universe (the density required for a flat 
universe, i.e & = 1). The density parameter, 2, essentially describes what the universe is 
made of by describing how large the contribution of a certain component is compared to the 
total energy density. Together these parameters for the different components of the universe 
decide the evolution and dynamics of the universe. Now using Equation 7 and the density 
parameter one can express Equation 5, and the Hubble parameter, as 


H(z)? = H3[O-(1+ z)* + On(1 + 2)? +a], (8) 


which is a useful expression when calculating cosmological distances (Barbara Ryden 2006b). 

In order to measure the density parameters one usually either observe how structure grow 
or by measuring distances. There are two main ways of measuring cosmological distances. 
One is using standard candles, which are objects with known intrinsic brightness, and you 
look at how the lightcurves have changed (for example redshifted, how much dimming etc). 
The other method is using object with known physical size, otherwise known as standard 
rulers. If one observes a standard ruler they can deduce the distance to them by looking 
at the angular size of the object in the sky (Barbara Ryden 2006b). In this paper we will 
mainly consider the standard candles. 


3.2. Measuring the expansion of the universe with supernovae 


One of the most useful type standard candles are Type Ia Supernovae (SNe Ia) due to their 
standardized lightcurves. A Type Ia supernova occurs when there is a close binary between 
a white dwarf and a companion star where matter is accreted to the white dwarf until it 
reaches Chandrasekhar limit of about 1.4 Mo (Mazzali et al. 2007). When this limit is 
reached the electron degeneracy pressure is not enough to overcome the gravitational pull 
and the white dwarf collapse. The increase in temperature and density results in an explosive 
fusion reaction which then destroys the white dwarf. The majority of the energy released 
is in the form of kinetic energy of the ejecta and is of the order of 1044J. The radiation 
energy, which is around 10*? J, arises from the fusion of °°Ni into radioactive cobalt and then 
eventually into stable iron (Karttunen H. et al. 2017). About 20 days after the explosion 
the SNe reaches it’s peak brightness, of about Mg = —19.3 magnitudes in the B-band 
(more about photometry in the next chapter), after which the supernovae quickly fades. 
Hamuy et al. 1995 and Phillips 1993 found that there is a correlation between the peak 
intrinsic brightness and the rate of decline following the maximum, namely that brighter 
supernovae have a slower decline. This along with the fact that brighter supernovae tend 
to be bluer (Tripp 1998) leads to the possibility of corrections of the lightcurve. There are 
several algorithms that performs these corrections to the shape and color of the lightcurves. 
Once this is done the SNe can be used as a standard candle and help constrain cosmological 
parameters along with other probes, such as the Cosmic Microwave Background and Baryon 
Acoustic Oscillations. 
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However there are several uncertainties when it comes to measurements of SNe and their 
use in cosmology. One of these uncertainties are the intergalactic dust the light from the 
SNe has to interact with on the way. 


3.3. Comoving frame of reference 


When discussing cosmological densities, such as the number density of galaxies, intergalactic 
dust grains, stars etc. it is advantageous to discuss them in a comoving frame of reference as 
will become apparent later. This is a frame of reference which expands with the expanding 
universe meaning that a particle which is at rest has a fixed comoving coordinate which 
is constant with time. Take this in contrast to a particle’s physical coordinates, which 
is the distance you would measure if you froze time and laid out a ruler. A particle’s 
physical distance increase with time due to the expansion of space. When converting between 
physical coordinates and comoving coordinates we are interested in the scale factor at the 
coordinates of interest. Thus using a comoving coordinate system is especially advantageous 
when analysing densities. In physical coordinates the density is proportional to the cube 
of the scale factor, a~?, meaning that through the history of the universe the expansion of 
space is likely the main contributor to the change of density. We can thus factor the change 
in density from cosmic expansion out by looking at a comoving coordinate system. In this 
system any change to the density is the result of some "real" astrophysical effect and not 
just due to the expanding universe. 


4 Photometry and photometric systems 


Because we are going to use data from SDSS which uses photometry and photometric filters 
to do their measurement we have to discuss photometry a bit. Photometric filters are used by 
telescopes in astronomical surveys to select the wavelengths to observe. The filter transmits 
one portion of the spectrum to the detector in the telescope and dismisses any other. In 
other words, they are band-pass filters. 

Transmission curves are a good representation of the behavior of these filters. Trans- 
mission curves essentially describe how efficiently the filter transmits light from a certain 
wavelength. There are two parameters that sufficiently describe the transmission curves. 
These are the effective wavelength (Acs fective) and the bandwidth, or otherwise known as 
the Full Width at Half Maximum (FWHM). Both of these parameters will be described in 
more detail in the following sections. 


4.1 Effective wavelength 


The effective wavelength is essentially the same as the isophotal wavelength (King 1952). 
The definition of isophotal wavelength is given in (Tokunaga and Vacca 2005) as: 

A source with the intrinsic spectral distribution F(A) cause N, number of photoelectrons 
per second inside the detector in a telescope. These photoelectrons are then detected and 
N, is described as: 

Np = ~ [rosea (9) 


where fh is Plancks constant, c is the speed of light in vacuum and S(A) is described by 
S(A) = TA)Q(A) RA) Arter (10) 


and is the total system response. Where T(A) describes the atmospheric transmission (i.e. 
disturbances from the atmosphere), Q(A) describes the quality of the telescope setup, R(A) 
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describes how good the filter is at transmitting light at a certain wavelength and A;,; is the 
area of the mirror collecting light in the telescope. Because of the mean value theorem for 
integration and the fact that Fy(A) and S(\) are both continuous and S(,) is always positive 
there exists a Ajgo such that 


Fy(Aieo) / MOIS 7 AFA(A)S(A)dA (11) 


where Ajso is the isophotal wavelength. If we then rearrange Equation 11 we get 


_ fAF\)S(A)dA 


FP) (Aiso) —_ f XAS(A)dr — (Fy) (12) 


where (F) is the mean value of the intrinsic flux of the source outside of the atmosphere 
across the bandwidth of the filter. This means that \;,. is the wavelength where the mean 
flux across the bandwidth equals the monochromatic flux F\(Aiso). This leads to Ais. being 
the wavelength that best represents the filter, or the transmission curve. Or in other words, 
“the instrument acts effectively as if all the light were concentrated at this wavelength” (quote 
from King 1952). As mentioned earlier js. is essentially the same as Ac f fective: 


4.2 Bandwidth (FWHM) 


The Full Width at Half Maximum (FWHM) describes the bandwidth over which the filter 
transmits the light. Meaning it is an interval of wavelengths where the lower and upper limit 
is defined as the wavelength where the filter achieves 50% of its maximum transmission. 
The light with a wavelength inside the bandwidth gets transmitted and the light with a 
wavelength outside gets dismissed. 


FWHM and Agr for the SDSS r filter 
0.5 





--- FWHM 
I ——- Aerrective 
! =—— SDSS r filter 
0.4 1 
! 
i] 
1 
I 
>0.3 4 
= 1 
> ! 
= | | Peapeee-see i a 
7) 1 
a 1 
H 0-2 : 
! 
1 
1 
0.1 : 
! 
i] 
I 
I 
0.0 ! 


5500 5750 6000 6250 6500 6750 7000 7250 
Wavelength [A] 


Figure 5: Illustration of the Full Width Half Maximum (FWHM) and the effective wave- 
length, Ae¢¢, for the SDSS r filter. 


4.3 SDSS filters 


There exists several different filters and photometric systems used in different surveys that 
want to look at different parts of the spectrum. Some may be interested in far infra-red 
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features and use appropriate filters for that, and some may be interested in UV features and 
hence use different filters. In this thesis we will use data obtained in the SDSS surveys. As 
such the data is given in magnitudes for the different SDSS filters. Below is a table containing 
the effective wavelength and bandwidth of the SDSS filters as well as a plot showing all of 
the transmission curves. 























Filter Aef fective [A] FWHM [A] 
u 3057 599 
g 4825 1379 
r 6261 1382 
i 7672 1535 
Zz 9097 1370 























Table 2: The effective wavelength and the FWHM for the different filters used in the SDSS 
survey as described in (M. Fukugita et al. 1996). 
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Figure 6: Illustration of the Full Width Half Maximum (FWHM) and the effective wave- 
length, Aer¢, for the SDSS r filter. 


4.4 Magnitudes 


Since humans perceive brightness, or change in brightness, on a logarithmic scale, astronomers 
have implemented the magnitude system. It is basically a logarithmic scale of the flux a star 
emits and it is calibrated to a certain star with a known brightness or some other calibration 
flux. Pogson’s definition of apparent magnitudes is defined by the equation 


m= —2.5logio(#), (13) 


where F' is the flux of the object and Fo is the calibration flux. Sometimes a calibration star 
is used (often Vega). Instead of Vega one can use several stars or some calibrated absolute 
unit of flux as in the AB magnitude system. When there are several filters looking at different 
wavelengths one often has to calibrate the magnitude separately for each filter due to the 
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difference in sensitivity for the different filters for the instrument. The magnitude system is 
designed such that a magnitude 1 star is exactly 100 times brighter than a magnitude 6 star. 
This means that the lower the magnitude the brighter the object is and each magnitude step 
changes the brightness by a factor of ¥/100 (Karttunen H. et al. 2017). One can then convert 
the fluxes observed through photometric filters to magnitudes for that corresponding filter. 
Notice however that the SDSS magnitudes differ slightly from Pogson magnitudes, especially 
for low signal to noise sources (M. Fukugita et al. 1996). 


4.5 Colour 


For any photometric system with multiple coloured filters one can define colour indices. The 
colour index is the difference in magnitude between two filters, taking the magnitude for 
the lower wavelength filter and subtracting the higher wavelength filter magnitude. Thus 
the smaller the colour index is the more blue the object is and conversely the larger the 
index is the redder the object is (Karttunen H. et al. 2017). Colours are a useful measure 
to compare observations to simulations with since it is a logarithmic difference of the fluxes 
of a spectrum. This allows removing any systematic factors which are applied evenly over 
the entire spectrum. Below is a derivation showing that any constant factor k that is evenly 
multiplied with the spectrum gets removed. 


Mz — My = 2.5logioFyk — 2.5logioFrk 
= 2.5(logioFy + logiok — logio Fx, — logiok) (14) 
= 2.5(logioFy = logioF x), 


where Equation 13 was used, x and y are two different filters and m is their magnitudes, 
Mz — My is their colour, k is the constant, and F' is their zero point corrected fluxes. 


5 Intergalactic dust and their influence on cosmological 
measurements 


The light emitted from distant objects may encounter several obstacles on the way, other 
galaxies, stars etc. However the majority of the time the light travels in the IGM, and 
on its path through the IGM the light may encounter dust grains. The existence of dust 
grains in the IGM has been a subject of great interest, and several studies have inferred a 
significant amount of intergalactic dust in the range of dust density Qaus: ~ 10~° — 107° 
(Loeb and Haiman 1997; Inoue and Kamaya 2004; Masataka Fukugita 2011). When it 
comes to intergalactic dust we are mainly interested in its interaction with light. There are 
several ways light can interact with dust. Dust grains can re-emit incoming UV and optical 
photons in the infrared to far-infrared. Dust grains are also able to scatter the X-rays from 
point sources with small scattering angles producing extended, diffuse halos instead of point 
sources. Additionally, some theoretical models predict the existence of small, rapidly rotating 
dust grains. If these dust grains have electric dipole moment there will be rotational emission 
emitting at microwave frequencies (Johansson 2015, B. T. Draine 2003a). 


5.1 Sources of dust 


The currently prevailing theories state that the producers behind dust are Core-Collapse 
(CC) supernovae (SNe) and Asymptotic Giant Branch (AGB) stars. As has been mentioned 
previously stars in their AGB phase undergoes strong stellar winds which can eject up to 
10-4Moyear~' in their late stages. During the AGB phase the star may lose up to ~ 80% 
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of their mass in the form of these stellar winds containing gas and dust. However it takes a 
long time for a star to evolve in to the ABG phase of its life, meaning that ABG stars are 
unlikely a contributor to dust at high redshift (Johansson 2015, Masataka Fukugita 2011, 
Gall, Hjorth, and Andersen 2011). 

Core Collapse SNe is a potential contributor to the existence of dust at high redshifts. 
CC SNe occurs for large stars which have a shorter lifespan thus CC SNe is not as limited 
by time as AGN stars. There have been several studies providing strong evidence that CC 
SNe produces dust, contrary to Type la SNe which has strong evidence not to contribute 
any dust production (Johansson 2015). 

To summarize: during the early period of the universe (towards high redshift) the main 
producers of dust is highly likely CC SNe, since the universe has not been around long 
enough for a significant amount of AGB stars to have formed. In contrast to later periods of 
the universe (lower redshifts) the main producers of dust will likely tip more towards AGB 
stars, since at later stages of the universe AGB stars become more abundant compared to 
massive stars that leads to CC SNe. This is partly because the lifespans of the two, where 
massive stars leading to CC SNe have a significantly lower lifespan compared to AGB stars. 
This results at lower redshifts the only CC SNe that occur are from relatively new born 
stars, while AGB stars are older stars coming from several different generations. Xiang et al. 
2018 looked into stellar mass distribution and star formation history of the galactic disc and 
confirms that there are more Mopc~® for older stars with ages over 1 billion compared to 
the ones younger than 1 billion. Another factor favoring AGB stars is the fact that stars 
tend to be less massive. It is very rare for stars to accumulate enough mass to incur a CC 
SNe (Karttunen H. et al. 2017). 


5.2 Intergalactic dust 


The existence of interstellar dust within galaxies has long been known. The existence of inter- 
galactic dust, however, is under more debate. There have been numerous studies conducted 
surrounding the existence of intergalactic dust and several of them have inferred roughly the 
same intergalactic dust density, Qgus¢ ~ 107° — 107°, from estimates of how stellar density 
and metallicity depends on redshift (Loeb and Haiman 1997, Masataka Fukugita 2011, Inoue 
and Kamaya 2004). 

There have also been studies surrounding the optical effects of intergalactic dust and 
observational proof of its existence, e.g reddening, extinction etc (short summary Johansson 
2015). 

One problem that intergalactic dust faces is that it must escape the potential well that is 
the gravitational field of the host galaxy where the originator star resides. There are several 
ways this could occur. The dust could be expelled from the originator with sufficient velocity 
due to radiation pressure such that they escape the potential well. During galaxy mergers 
interstellar dust may be expelled into the IGM or due to ram pressure as the galaxy moves 
through the IGM (Bianchi and Ferrara 2005; Aguirre 1999). This results in the composition 
of intergalactic dust differ from that of interstellar dust since intergalactic dust are more 
weighted towards heavier grains. Larger grains often result in less reddening and higher 
values of R,, and are often labeled as "grey" dust. 


5.3 Potential dust types 


From examining meteorites and investigating how dust grains interact with light, one can 
deduce what these grains are made of. Draine (B. T. Draine 2003a) made a short list of 
likely candidates for interstellar dust grains based on the observations just mentioned. It was 
concluded that silicates, which are any anion consisting of silicone and oxygen, are a major 
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contributor to the total mass of interstellar dust. Carbonaceous materials, or materials which 
mainly consist of carbon, were also a prominent group of materials, in particular graphite as 
indicated by the 2175 A feature in the spectrum of quasars. There are also other materials 
like SiC and carbonates which is not as prominent (B. T. Draine 2003a). J. S. Mathis, 
Rumpl, and Nordsieck 1977 found that graphite was a necessary component to any dust 
grain mixture, in order to achieve the correct spectral features. Graphite along with any 
other silicates proved to be a successful mixture in reproducing the interstellar extinction 
of the Milky Way, with a size distribution of graphite from 0.005 wm to 1 wm and silicates 
from 0.0254m to 0.25 um. This distribution followed the power law dn/da « a~* with 
k ranging between 3.3 — 3.6. However since intergalactic dust is weighed towards larger 
grains, as discussed in the previous section, we find a different dust size range for those. A 
intergalactic dust model that was found to be viable had a distribution of either silicate or 
graphite, or a mixture of both, with grain sizes in the range 0.25 zm to 0.05 xm. So for our 
model we assume a composition of either silicate or graphite. 


5.4 X-ray scattering 


It has been shown that intergalactic dust grains may produce halos of scattered light around 
X-ray sources on cosmological distances due to X-ray scattering. The properties of the halo 
may provide us with information about the intergalactic dust causing the halo (John S. 
Mathis and Lee 1991, Smith and Dwek 1998, B. T. Draine 2003b). The differential cross 
section of the scattering interaction can be calculated using the Rayleigh-Gans approximation 
for a spherical dust grain with the radius a using the following equation (from Hayakawa 
1970; Mauche Gorenstein 1986): 








Cece = Ap (Se) [AE] tense) (15 


Here x« = (4ra/X)sin(0/2), 71 is the spherical Bessel function of the first order j;(”) = 
(sina) /x? — cos(x)/x and Ag is a normalization constant which depends on the energy E of 


the X-rays interacting with the dust grains and the grain composition. Ag can be calculated 
using the following equation: 


v= 1a(Gr) (gems) (AE ) a 





where Pgrain is the mass density of the dust grains, Z and M is the atomic number/charge 
and mass number respectively and FF) is the atomic scattering factor which depends on 
the X-ray energy (Henke, Gullikson, and Davis 1993). We see that Equation 15 only depends 
on the wavelength/energy of the X-rays, the scattering angle, and dust grain properties. If 
we know the source of the X-rays and their energies and look at the halo produced to deduce 
the scattering angle we are able to gather some information about the properties and amount 
of the dust grains which the X-rays interact with. The central core of this equation is also 
approximately Gaussian with the root mean square (rms) of the scattering angle indicating 
typical sizes of the scattered halo, 0-ms * 62.4(1.0um/a)(1keV/E)arcsec (Johansson 2015). 


5.5 Dust extinction 


Light might undergo extinction due to scattering or absorption of the light when it encounters 
dust grains (see Figure 10). Light which has been affected by extinction are often linearly 
polarized since extinction are polarization-dependant requiring that some dust grains are 
aligned with the local magnetic field. Often extinction leads to reddening of the light, 
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which is an important factor when it comes to standard candles, where redshift is used as a 
distance measurement. This means that extinction of the light produces an offset of distance 
measurements towards greater distances. Thus, in order to make accurate measurements one 
has to take into account light extinction due to dust in the line of sight (Johansson 2015, 
B. T. Draine 2003a, Johansson and Mortsell 2012). 

The total to selective extinction ratio can be estimated by the ratio: 


A Ay 


Ry = = 
v=" E(B-V) Ap-— Ay 





(17) 


where £:(B—TV) is the colour excess, and Ay and Ag is the extinction in the V- and B-band 
respectively. Ry could be seen as the "reddening parameter" or the slope of the extinction 
curve in the optical region and can be used to examine how much extinction has occurred 
(Johansson 2015, Johansson and Mortsell 2012). When examining the effect of dust on 
measurements through dust extinction there are two key factors to take into consideration. 
One is the number density of dust grains between the observer and the object of interest. 
The other one is the extinction cross section, e.g with a given amount of dust grains how 
likely is it that a photon will interact with a dust grain on its path. 


5.5.1 Comoving number density of dust 


Since intergalactic dust spawns from stars we make the assumption that we can derive the 
comoving number density of intergalactic dust from the comoving star formation rate. This is 
one of our main assumptions of the model. The star formation rate, 6,(z), given in Hernquist 
and Springel 2003 is used. We then introduce a normalization factor, no, such that 


= Ot ___Px(z’) 
n(z) = no | dz (it 2e(z)’ (18) 





where ¢(z) = ae = /0,(1+ 2)4 + On(1 + 2)? + Qa is the dimensionless Hubble param- 
eter from Equation 8 (note that ©, is very small in comparison to Q,, and Qa, making 
it negligible). We essentially normalize the integral such that no is the comoving number 
density at redshift 0. We can thus relate no to the present dust density of the universe Ogust 


through (Johansson and Moértsell 2012) 
Quaust = (no *2a* pgrain)/ Perit, (19) 


where a is the radius of the dust grain which is assumed to be spherical. pgrain is the density 
of the dust grain and pepiz = 1.88x10~79h? g cm~® is the critical density of the universe. We 
can now rearrange the above equation to 


2 Qudust 3Perit 


No m 
a? AT Pgrain 
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Now for a given Qays_ we can calculate the no for a particular dust size and dust composition. 
In Figure 7 we see the normalized comoving number density of dust. 
5.5.2 Dust scattering & absorption cross section 


The geometrical cross section is the probability of a photon interacting with a dust grain, 
taking into account the size of the grain and the probability of dust-photon interactions. The 
cross section for a particular wavelength, A, can be calculated as 


o(A) = ra?Q(A), (21) 
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Figure 7: The evolution of the comoving number density of intergalactic dust. As can be seen 
there is a greater amount of dust towards lower redshifts. This should come as no surprise 
since the universe is continuously enriched with dust. 


where a is the radius of the dust grains which are assumed to be spherical, Q is the sum of 
the efficiencies of the interactions due to scattering and absorption. Throughout this thesis 
a composition of silicates and graphite will be assumed for the intergalactic dust. From J. S. 
Mathis, Rumpl, and Nordsieck 1977 we know the sizes of dust grains, and the sum of the 
efficiencies from Laor and Bruce T. Draine 1993 is used for silicates and graphite. In Figure 
8 we see the data from Laor and Bruce T. Draine 1993 showing the silicate and graphite 
efficiencies. Not that these are laboratory results and not astronomical measurements. 


5.5.3 Optical depth of intergalactic dust 


Given the previous discussion we can now estimate the optical depth. The total optical depth 
due to scattering and/or absorption due to interactions with dust grains between redshift 0 
and Zp; can be calculated with the following formula: 


cpm neo + 2)? , 
al a(z") dz’, (22) 


where o(X’) is the geometrical cross section and n(z’) is the comoving number densities of 
the dust grains from the previous sections and \’ = A,ps/(1+ z). For the sake of this thesis 
we assume dust to consist of homogeneous sizes, i.e constant size for each model. 

The optical depth of a medium is the logarithm of the ratio between observed to trans- 
mitted intensity of light as it has traveled through that medium. It essentially describes how 
much light passes through a certain medium, and how much light is lost on the way. How 
the optical depth affects the observed intensity is described by 








T(Aobs; Zem) = 


Iops(A) = Ipe~", (23) 


where I,y, is the observed intensity and Jp is the intrinsic intensity of the observed object. 
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Figure 8: Extinction efficiencies (Qeztinction = Qscattering + Qabsorption) for both silicate and 
graphite for varying dust grain sizes over wavelength. As well as the cross section for both 
grain materials. Data for the efficiencies are from Laor and Bruce T. Draine 1993 


With the optical depth one can also calculate the observed dust attenuated flux F%,, of 
a certain object at a redshift z-;, observed at a wavelength A.», using 


Fops(Aobs: Zem) = Pom : eW Text (Aobe,Zem), (24) 


where F.,, is the intrinsic flux, Tez; is the optical depth due to extinction which can be 
calculated using Equation 22. This means that the dimming of light emitted at z.,, and 
observed at wavelength A,,, can be calculated as (Johansson 2015, Johansson and Mortsell 
2012) 


2.5 


Am(Xobs; Zem) = Trio reat obs Ber). (25) 


6 Using quasars as back light sources 


In the previous chapter we discussed the light transmission properties of intergalactic dust. 
In order to measure these properties from observations we require a well determined back 
light source enabling us to detect the transmission of its light as it passes through the 
intergalactic dust. As we will see in the following description quasars are particularly well 
suited astronomical sources for this task. 
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QUASAR ENGINES 


The accretion of matter onto giant black 
holes in the centres of galaxies powers 
the extreme luminosities of quasars. 
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Figure 9: An illustration of an AGN, showing the supermassive black hole in the center with 
the accretion disc around it. The dust torus that surrounds the black hole is also shown as 
well as the particle jets. It also shows how the viewing angle change our definition of the 
object. Image source: Antonucci 2013 


6.1 A Quasi-Stellar Object (QSO) 


The current state of knowledge says that quasar, or Quasi-Stellar Objects (QSO’s), is an 
extremely luminous type of Active Galactic Nuclei (AGN). AGN’s are compact regions in the 
centers of galaxies with greater than normal luminosities with characteristics indicating that 
the luminosity is not coming from stars. Sometimes the AGN may be brighter than the entire 
galaxy it resides in. The leading models suggest that supermassive black holes with masses 
between 107 — 10!° Mg reside in the centers of galaxies (Cattaneo et al. 2009; Karttunen H. 
et al. 2017). Around the black holes there is an accretion disk of gas and dust, and beyond 
the accretion disk lies a massive dust torus supplying the accretion disc with materials. It 
is believed that the source of the energy for an AGN comes from the gravitational potential 
energy of this accretion disc as the gas and dust falls into the black hole. In some cases this 
may produce jets ejecting matter up to relativistic speeds perpendicular to the rotational 
plane of the accretion disc. See Figure 9 for an illustration of an AGN. 

The conversion of gravitational potential energy into radiation is a more efficient way of 
converting mass into energy compared to regular fusion in stars. The efficiency ranges from 
5.7% for a non-rotating Schwarzschild black hole to 32% for a rapidly rotating Kerr black 
hole. Comparing this to the efficiency of nuclear fusion of hydrogen of about 0.7% we see 
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that gravitational energy released by accretion is far superior in its efficiency (Lambourne 
2010). Due to the large energy production of AGN galaxies it is unlikely that it will continue 
through the entire lifespan of the galaxy, and the AGN will eventually quiet down as the 
mass in the accretion disc is depleted. Thus AGN’s are likely not a separate type of galaxies 
but just a phase that only some galaxies goes through before settling as "regular" galaxies 
(Karttunen H. et al. 2017). 

The high luminosity of quasars are believed to be due to the viewing angle of the AGN 
(Karttunen H. et al. 2017). As seen in Figure 9 if viewed close to the particle jet the AGN 
will appear as a quasar, in contrast to when viewed from the side through the dust torus 
when the AGN appears as a radio galaxy. If the particle jet were pointed directly, or almost 
directly at us the AGN would appear as a strong radio source due to the particles in the 
jets move at relativistic speeds. These objects are known as a blazars (Karttunen H. et al. 
2017). 

It has been noted that quasars vary greatly in luminosity with time. Sometimes on the 
range of a couple of days. This means that the quasars generate and emit the light from a 
small region, since the entire region has to be causally linked on those timescales to incur 
a coordinated luminosity change. So when observing quasars they appear almost like point 
sources, which is why they originally were mistaken to be a type of star (Karttunen H. et 
al. 2017). 

Over 750 000 quasars have been observed as of August 2020 (Lyke et al. 2020), most of 
which has been found by SDSS. It has been shown that quasars were more common in the 
early universe, i.e towards larger redshifts (Karttunen H. et al. 2017). This might be due 
to the energy production of quasars not being sustainable on large timescales as well as the 
universe being more active in the past (see Section 2). 


6.2 Properties of QSO’s 


Quasars are detectable over the entire observable electromagnetic spectrum, from radio to 
X- and gamma rays. Generally quasars are brightest around their rest-frame ultraviolet 
wavelength of 1216 Angstréms, which is the Lyman-alpha hydrogen emission line. However 
many of the observed quasars are at larger redshifts (Lyke et al. 2020) meaning that in the 
frame of the observer the peak luminosity is often towards larger wavelengths, sometimes 
even towards the near infrared part of the spectrum (Karttunen H. et al. 2017). 

For many quasars their spectra can be described as a featureless continuum in the ultra- 
violet to optical region with some strong and broad emission lines. When comparing several 
quasars to each other it has been apparent that most quasars are strikingly similar in terms 
of their spectra. It has been shown that quasars are similar in terms of both spectral features 
as well as colours. These similarities even exist over large ranges of magnitudes in their lu- 
minosities as well as ranging redshift (Vanden Berk et al. 2001; Telfer et al. 2002; Johansson 
and Mortsell 2012). To analyze the difference in emission properties, i.e spectra, studies have 
looked at differences in the rest equivalent widths of some spectral lines such as Lya, as well 
as differences in spectral slope (Peng et al. 1999). The spectral slope is essentially the angle 
of the slope of the spectrum in a logarithmic coordinate system. 

To the extent that quasars are homogeneous in terms of spectra is found in no other 
astronomical object, which is why quasars are a great object to use as a back light source. 
Due to the relative homogeneity of quasars one can construct composite template spectra of 
quasars from a sample of real quasar spectra, i.e an average spectrum (Vanden Berk et al. 
2001) (more on this in Section 7.2). Using this composite spectrum one can study general 
properties of a quasar population, which in our case is the effect of intergalactic dust. 

To do that we include this template spectrum in different dust models, where we apply 
our dust attenuation model to the composite spectrum and compare the results to obser- 
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vations. We are then able to compare the difference between the result of the model and 
the real observations for each model. For this thesis we are going to compare the models to 
observations using colours (more on this in Section 7.5). Due to the homogeneity of quasar 
colours the difference between the models to observations will be due to differences in the 
models and not because of differences in the quasars. 

To conclude, here are some of the main points why quasars are a good astronomical object 
to use as our back light source: 


e Quasars are relatively homogeneous in terms of both spectral features and colours 
(Johansson and Mortsell 2012). 


e There are a lot of measurements of quasars over a large span of redshifts (Lyke et al. 
2020). 


e Quasars are comparable over the entire redshift span (Vanden Berk et al. 2001). 


e There exist composite spectra which could be used as the simulated back light source 
(Vanden Berk et al. 2001; Telfer et al. 2002). 


7 Method & Result 


This master thesis builds on previous studies of the effects of intergalactic dust. Much of the 
method is inspired by Johansson and Mértsell 2012; Mértsell and Goobar 2003; and Ostman 
and Mortsell 2005. Below are descriptions of the methods used, thoughts behind each step 
and the resulting figures. The cosmological parameters of the Planck Collaboration et al. 
2020 is used. 


7.1 The experimental setup 


In this thesis we are interested in the influence of intergalactic dust on observations. To 
investigate this we construct a forward model of the effects of intergalactic dust on light 
emitted by quasars following the conceptual outline shown in Figure 10. We construct the 
model such that we have a simulated distant quasar that emits light into the intergalactic 
medium where the photons interacts with intergalactic dust grains getting attenuated by dust 
through scattering and absorption (see Figure 10). How the emitted light interacts with the 
intergalactic dust is governed by the known equations explained in Section 5.5. We have now 
modelled our transmitted quasar light and can now compare our model to observations. This 
comparison will be done using colours, which will later be shown to be optimal for this thesis. 
If the simulations and real observations agree in terms of the results it means the model is 
a good representation of real intergalactic dust. And if they do not agree we can tweak the 
model, like changing the composition or the size of the dust grain sizes. In doing so we are 
able to constrain the model to fit real observations thus calculating the most probable dust 
scenarios. In order to do all of this we must first find a way to simulate the quasars as our 
back light source. 


7.2 Simulating quasars as the back light sources 


As we saw in Section 6 quasars are a great source to use as our back light source. An 
additional advantage is that there exist photometric observations for a lot of quasars even 
towards higher redshifts (z ~ 3) which is advantageous since we also want to investigate the 
redshift dependence of our dust model. Now in order for us to be able to apply our dust 
attenuation model onto our back light source we need a simulated spectrum for our back 
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Figure 10: An illustration of the experimental setup of the model. A quasar emitting light 
through intergalactic dust. As the light travels through the dust it undergoes both scattering 
and absorption, causing a reddening and dimming of the emitted light as it reaches the 
observer. 


light object. For the simulated spectrum we can use a composite spectrum of quasars. We 
use the median composite spectrum derived in Vanden Berk et al. 2001, which used over 
2200 samples of spectra from SDSS. This composite spectrum spans from 800 to 8555 A. We 
splice this spectrum with the composite spectrum produced in Telfer et al. 2002 which used 
the spectra from 101 quasars observed with the Hubble Space Telescope. This spectrum 
covered the wavelengths 309 to 3386 A. The splicing of the spectrum was done to achieve 
a greater wavelength span of the resulting composite spectrum. We disregarded the part of 
the spectrum below 342 A due to the noise and larger errors caused by the Lyman alpha 
forest. We spliced the spectra at about wavelength 1100 A to produce the total composite 
spectra. The resulting composite spectrum can be seen in Figure 11. 


7.2.1 The composite spectrum 


As previously mentioned in this thesis we will use composite spectrum as our simulated back 
light source spectrum. A composite spectrum is the mean spectrum from several observed 
spectra, including as many spectral features of the source as possible. There are a couple of 
steps to creating a composite spectrum (Vanden Berk et al. 2001): 


1. Selecting the input spectra used to generate the composite. 
2. Determine accurate redshifts of the spectra. 

3. Shifting the spectra to our rest frame. 

4. Scaling or normalizing the spectra. 
5 


. Divide the spectra into wavelength bins and calculate the mean spectra for each bin 
then combine them to generate the complete composite spectrum. 


To do this the SDSS composite spectrum (Vanden Berk et al. 2001) used SDSS spectra data 
for 2204 quasar spanning the redshift range 0.044 < z < 4.789 with a median quasar redshift 
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Figure 11: The complete quasar composite spectrum in arbitrary flux units. 


of z = 1.253. The flux calibration for these spectra used several F subdwarf stars to define 
the response function. This largely corrects for Galactic extinction. Other steps were also 
made in the data selection process to minimize the effect of host galaxy extinction (Vanden 
Berk et al. 2001). 

Vanden Berk et al. 2001 created two different composites, the first which preserves the 
spectral features of the combined quasars, and the second which preserves the global shape 
of the continuum of the combined quasar spectra. For this thesis we use the first composite 
which preserves the relative fluxes of the emission features. This is because we are going to 
be looking at quasar colours, thus the relative fluxes depending on wavelength are important. 

In order to determine accurate redshifts Vanden Berk et al. 2001 created a smaller com- 
posite spectrum focusing on the [OIII] emission line at wavelength 45007 A which could be 
used as a template. They constructed the composite using 373 QSO spectra with strong 
[OIT]] emissions lines. The wavelength array for each spectra was then shifted to the rest 
frame based on the position of the [OIII] line and the spectra was rebinned into bins of 1 
A width in the wavelength space while preserving the flux. The flux density of the bins 
was then scaled and normalized depending on redshift (Vanden Berk et al. 2001). The fi- 
nal [OIII] composite spectrum was then generated by summing the flux density of each bin 
for the shifted, rebinned, and scaled spectra. Using the [OII]] composite they were able to 
more accurately define the redshifts of the spectra missing a strong [OIII]-line or which had 
the [OIII]-line shifted beyond 9200 A. Then the process behind the creation of the [OIII| 
composite was done again for the complete composite spectrum using all of the 2204 QSO 
spectra. To create the median composite spectrum preserving relative fluxes they took the 
median flux density for each bin from the shifted, rebinned and scaled spectra. 

The HST composite spectrum (Telfer et al. 2002) is generated in a similar fashion as 
explained above but using a different data set focused more towards bluer wavelengths. 
They used 322 HST spectra of 184 QSO’s with z > 0.33 in order to study ultraviolet spectral 
properties of QSO’s. 

A problem with using these composite spectra for our analysis is that of course the spectra 
used in the generation of the composites might also be attenuated by intergalactic dust. 
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Figure 12: A flowchart showing the iterative correction work that was done to correct for 
potential intergalactic dust effects in the composite spectrum done in Mértsell and Goobar 
2003 


7.2.2 Possible contamination in the generation of the composites 


Much has been done in the creation of the composite spectrum to account for Galactic dust 
extinction and dust extinction from the host galaxy. However there could be cases where the 
observed spectra they used to generate the composite already are attenuated by intergalactic 
dust. This is clearly a nuisance since we would prefer the spectrum to be free of dust effects 
since we are applying a dust attenuation model onto the spectrum. 

Mortsell and Goobar 2003 formed an iterative correction procedure to investigate the 
influence of possible contamination. They began by using the composite spectrum from 
Vanden Berk et al. 2001 and using the known dust attenuation equations to guess the dust 
contribution. Given this guess they were able to correct the original spectrum for the guessed 
dust contributions. With the new composite spectrum they were able to redo the process of 
guessing the dust contribution then they iteratively did this entire process until the composite 
spectrum converged. This process thus reduces the potential effects of intergalactic dust. A 
flowchart showing this iterative process can be seen in Figure 12. 

Now using the final dust-corrected composite spectrum they did the initial analysis again 
and compared the result with the results from the non-corrected composite. They found that 
the possible effects from intergalactic dust did not have a sufficiently significant effect on the 
result in their 7 analysis for the result to be invalidated. However this is a systematic effect 
that is important to take into account. Since their analysis is very similar to ours where we 
both used quasar colours to analyse intergalactic dust properties with a x?-analysis we also 
assume these effects would have little influence on our results. 

Studying the potential systematic effects in the composite spectra is of great importance 
in improving the statistics of the result. As will be seen in Section 7.8 the error in the results 
is governed by the errors in the composite spectrum. Reducing the systematic effects as well 
as improving the statistics of the quasar surveys should improve the control of these effects 
when generating the composites and thus improve on the error of the composite. 
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7.3 Applying dust attenuation to the back light sources 


Now we want to apply the dust attenuation, caused by the photons travelling through the 
intergalactic dust, to the spectrum of the source. For a certain dust model we are now able to 
firstly redshift the spectra to the redshift we want to look at, then apply the dust attenuation 
to the flux of the back light source spectra. This is done at each wavelength of the spectrum 
using Equations 22 and 24. The cross section in the optical depth integral was easily attained 
from the data from Laor and Bruce T. Draine 1993 for both silicate and graphite, i.e the 
materials we are assuming our intergalactic dust to consist of. For the comoving number 
density of dust we used two different models, or approaches. 


7.3.1 Constraining the total amount of all dust grains 


Initially we just related the comoving number density of dust to the comoving star formation 
rate, without the inclusion of the no factor from Qaus¢ (see Section 3.3). 


ey _ bxl2’) 
n(z) « | dz (i+ ze(2): (26) 
This is an additional simplification of the model, essentially constraining the total amount of 
dust grains in the universe to a constant. This is because we have no grain size dependency 
in this number density, thus changing the grain size of this model does not change to total 
amount of dust grains. From now on this will be referred to as Model 1. Model 1 was the 
first model tested since it is an additional simplification of the model, making the model 
less complex since we have no constrain from OQgust. It was also made to investigate certain 
dependencies, for example how would the effect of dust change if all dust grains increased /de- 
creased in size but not change in amount. This could help with checking the validity of the 
algorithm, seeing if certain steps produce the expected result etc. 

Since for this model the comoving number density is invariant under dust size changes 
it results in tT « a? due to the only dust size dependency in the optical depth is from the 
cross section which has a a? dependency. This means when looking at how the emitted flux 
is altered by the dust attenuation (See Equation 24) we see that Fou, x e-®. This is quite 
an unstable dependency on the observed flux since small changes in the dust grain size will 
cause large changes in the observed flux. 

All of the plots produced for Model 1 in the following sections can be seen in Appendix 
A. 


7.3.2  Constraining total mass of all dust grains 


The second model included the constrain on the total amount of dust grains from OQdgust by 
implementing the normalization factor no from Equation 20. Instead of constraining the total 
amount of dust grain in the universe to a constant this model constrains the total combined 
mass of all dust grains in the universe. Thus if we increase the size of the dust grains we 
reduce the total amount of dust grains. This is because of the new dust grain size dependency 
in the comoving number density, n « a~?, that comes from the constrain from Qgys¢. From 
now on this will be referred to as Model 28. Model 2 is a more physically appropriate 
model since putting constraints on the total combined mass of dust grains makes more sense 
from a scientific perspective as compared to putting constraints on the total amount of dust 
grains. All of the plots shown in the method and results section will be using Model 2 with 
Qaust = 5 X 10~® (as given in Ménard et al. 2010). 

Since this model has the n « a~? dependence in the comoving number density we get the 
optical depth dependency of t x a~! since now we have both the dust grain dependency in 
the cross section, a?, and in the comoving number density, a~°, which are multiplied with 
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each other. This leads to the observed flux being proportional to the dust grain size through 
Fyps x e7/ “, which is not as sensitive to changes in the dust grain compared to Model 1. 
A summary of the features of the two different models is seen in Table 3. 
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Table 3: A summary of some of the main features and differences of the two models used in 
this thesis. 


7.3.3 Integrating the line-of-sight optical depth of intergalactic dust 


A major part of this thesis has been about making the calculations efficient and not too time 
consuming. This was especially important when it comes to calculating the optical depth 
integral (Equation 22) since, as will be discussed later, this integral will be calculated many 
times over. Even calculating the integral once requires a lot of calculations. This is partly 
because the comoving number density factor in the integral is itself an integral (Equation 
18). This means that at every point during the calculation of the optical depth integral we 
have to calculate a comoving number density integral. 

Because we are going to do a lot of calculations using the optical depth integral the 
calculation time had to be shortened. One way to do this is to interpolate the parts of the 
integral. Interpolations was used for the comoving dust density and the cross-section. An 
issue with using interpolations is that it reduces the accuracy of the calculations somewhat. 
This depends on how many points we are interpolating. For the cross-section 50 points were 
used and for the comoving dust density 200 points were used. The reason why we could 
take 200 points for the dust density will become apparent later on. We still felt that the 
increase in error was worth the reduced computation time since the computation time got 
reduced by a factor of 100 and the numerical error from the integral calculation went from 
a 0.000113% to 0.0024%. Note that this is not the total error, but only the error from the 
integral calculation. In Figure 13 we see the result of the optical depth integral ranging over 
redshifts in the range 0 < z < 3. The first plot is for varying dust grain sizes, and the second 
is for varying wavelengths. Both of the plots are for silicates. 


7.3.4 Predicting dust attenuated quasar spectrum 


According to the forward model described in Section 7.1, we will use the composite spectrum 
to estimate the emitted flux of a quasar at redshift z. To do this we first redshift the spectrum 
to the quasar’s restframe. This is done by multiplying the wavelengths by 1+ z then at each 
redshifted wavelength of the composite spectrum we multiply the flux with e~Tet Qevs:%em) 
where Text(Aobs; Zem) is from Equation 22. Practically this is done using interpolations. The 
optical depth is calculated for the redshifted wavelength array of length 50, equally spaced, 
covering the entire span of the redshifted spectrum. Since the optical depth’s dependence on 
wavelength is rather smooth using interpolation does not incur that great of an error. This 
means that for each spectrum we want to attenuate we have to calculate the optical depth 
integral 50 times, once for each wavelength bin. 

We see how the optical depth changes over wavelength for both varying dust sizes and 
redshifts for silicates in Figure 14. Now that we have Tert(Aods, Zem) we can then apply 
the dust attenuation on to the entire composite spectrum to produce our simulated dust 
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Figure 13: The optical depth as a function of redshift, r(z). First plot with constant wave- 
length and varying dust size and the second plot with constant dust size and varying wave- 


length. Here we note increasing optical depth with increasing redshift, as well as with 
decreasing dust grain size. 


. Optical depth's A dependance with varying redshift Optical depth's A dependance with varying aqgust 
10' 


Text(Aobs) 
Text(Aobs) 





104 


Wavelength [A] Wavelength [A] 


Constant dust size, a = 0.01 wm Constant redshift, z = 2 


Figure 14: The optical depth as a function of wavelength, 7(A.»;) for constant dust grain size 
(left figure) and with constant redshift (right figure). As seen in the left figure with increasing 
redshift the optical depth increase as expected since the light has to travel through more dust. 
The the shape of the curves in the right figure is caused by the extinction efficiency which 


differ for different dust grain sizes. Generally we have larger dust attenuation towards lower 
wavelengths. 


attenuated spectrum. In Figure 15 we can see a couple of spectra which has been affected 
by dust attenuation with different dust sizes or redshifts for silicate dust grains. 
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Figure 15: The predicted intergalactic dust attenuated spectra of back light sources with a 
constant dust grain size of a = 0.01 um with changing redshift (left figure) and with constant 
redshift z = 2 with changing dust grain size (right figure). Both are for silicates. In the left 
figure we see the spectra are shifted in the wavelength span due to their redshifts, as well as 
the fluxes being attenuated more towards larger redshifts at lower wavelengths. In the right 
figure we see larger dust attenuation for smaller dust grain sizes towards lower wavelengths. 
This is caused by the dependencies shown in Figure 14 


7.4 Discussion of the dust attenuation model 


Now lets discuss the the application of the intergalactic dust attenuation through the optical 
depth integral. We start by looking at the comoving number density. 


7.4.1 The distribution and number density of intergalactic dust 


For the number density we used two different models (see Sections 7.3.1 and 7.3.2). We choose 
to use Model 2 throughout this thesis (except for Section 7.6) since it is more plausible for 
the number density of dust to be correlated to the total amount of dust. 

Model 2 has the constrain that with changing dust grain size we set a constant mass 
density of dust grains in the universe, Qays:. This creates a relation between the dust grain 
size and the amount of dust grains. Increasing the dust size decreases the amount of dust 
grains. Qgust is something we can constrain by observations, while if we look at Model 1 
where with changing dust grain size we constrain the total amount of dust grains to be a 
constant. Here as we increase the dust grain size the total amount of dust grains will remain 
the same, resulting in the amount of dust in mass units increase. Essentially here we have a 
constant total amount of dust grains, which would make little sense it if were independent 
on dust grain size. It is also more difficult to constrain the amount dust grains, independent 
on size, through observations . Since what actually causes the effects from dust depends on 
the mass density of the dust grains, and not strictly the amount. 

For both of the models we assume that intergalactic dust is produced by stars and expelled 
from their galaxies. To do this we assume that the number density of dust follows the star 
formation rate given in Hernquist and Springel 2003. However this is a slight simplification 
since not all of the dust is created during the formation of stars. Instead they are created 
during later stages of their evolution (see Section 2. Another simplification we make for 
the number density of dust is that we assume the distribution of intergalactic dust to be 
isotropic, i.e the same in all directions, which is not strictly true. Since intergalactic dust 
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Figure 16: A conceptual image showing the possibilities of different quasars being more or less 
influenced by intergalactic dust compared to others. These differences depends on the spatial 
distribution of intergalactic dust. For the quasar emitting light through intergalactic dust we 
see the light gets both reddened and dimmed, compared to the light from the un-attenuated 
quasar. 


spawns from galaxies we expect higher concentrations of dust around galaxies and galaxies 
are not evenly distributed in space (Barbara Ryden 2006a). We would thus expect higher 
number densities of dust towards regions with more galaxies compared to regions with less 
galaxies. This means that not all quasar spectra at equal distance are affected equally by 
intergalactic dust. Some quasar light may travel through denser dust regions compared to 
others, as depicted in Figure 16. One could take advantage of this fact when generating 
composite spectra to reduce the systematic errors from by possible contaminations from 
intergalactic dust. To do this one could generate a intergalactic dust distribution map from 
the distribution of galaxies, then choose spectra from quasars which are in line of sights with 
relatively low amounts of dust to generate the composite. Another way to produce a cleaner 
composite spectrum would be to create a composite spectrum from quasar observations at 
lower redshifts. Since the further away the object is the more intergalactic dust the light will 
encounter on its way to us. By reducing the distance of light travel it would also reduce the 
effects from intergalactic dust and thus reduce the systematic errors (Ostman and Mértsell 
2005). 

Taking into account the spatial distribution of dust one could also produce more accurate 
dust models since we would have more accurate numberdensities of dust when we compare 
the models to real observations. We could also use this to more accurately correct for the 
effects intergalactic dust have on observations when taking into account how much dust is 
in the line of sight of each observation. The inhomogeneity of the dust distribution was 
something this thesis originally aimed at analysing, however other parts of the thesis took 
over and creating the algorithms took longer than expected. 

Assuming a homogeneous distribution of dust likely does not incur that great of an 
error, since we essentially take the average of all directions and use observations from all 
directions. However with increasing precision of measurements such small systematic errors, 
as assuming a homogeneous dust distribution, may become of larger importance. We see in 
Figure 7 the comoving number density of dust using the star formation rate in Hernquist 
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and Springel 2003. Here we notice the increase in number density towards lower redshifts, 
which is expected since the intergalactic medium is continuously enriched by dust through 
various processes. 


7.4.2 The properties of dust, their structure and cross section 


When it comes to the actual intergalactic dust grains several simplifications were made. 
Firstly for our model we only assume either silicate or graphite, while in reality it is likely a 
mixture of them both with some smaller contributions from other materials (B. T. Draine 
2003a). We also assume only spherical dust grains with the same dust grain size for the 
entire dust population. The data used for the interaction efficiencies are also from laboratory 
measurements and not astronomical measurements (Laor and Bruce T. Draine 1993). In the 
end I believe the assumptions made are valid and the same assumptions are made in several 
previous studies. It has been shown that a composition of spherical grains of either silicate 
or graphite make up for a valid dust model (Johansson and Mortsell 2012). 


7.4.3 The dependency of optical depth on model parameters 


Now we come to the result of the optical depth integral. As seen in figure 13 we see the 
evolution of the optical depth over redshift for varying grain sizes and wavelengths. We see 
the general trend of the optical depth increasing with redshift, which is expected. We also 
see that the optical depth increase with decreasing grain size, as expected considering the 
7 xa‘ dependency of the optical depth for Model 2 (see Section 7.3.2). 

We can now take a look at how the optical depth evolves over wavelength which is 
important since it determines how the optical depth is applied to the composite spectrum (see 
Section 7.3.4 which we will discuss in the next section. In Figure 14 we can see the evolution 
of optical depth over wavelength. In the left figure we see the same redshift dependency 
as discussed above, with increasing redshift the optical depth increases. In the right figure 
we see that the evolution of optical depth is highly dependant on the dust grain size. This 
is because the only wavelength dependency we have in the optical depth comes from the 
extinction efficiency, which is also highly dependant on dust grain size as seen in Figure 8. 
The shape of the optical depth curve in these plots is thus caused by shape of the extinction 
efficiencies. This is further supported by the tests done in Section 7.6 where we concluded that 
the only wavelength dependency in the optical depth comes from the extinction efficiency. 


7.4.4 Dust attenuated quasar spectra 


We can now discuss the predicted dust attenuated spectra, where we applied our dust atten- 
uation model onto the composite spectra. Here we essentially just multiplied the redshifted 
composite spectra with e7**(#) where Tex+(z, A) is the functions shown in Figure 14. We 
see the result of this in Figure 15, where we see with varying redshifts firstly we see the shift 
of the spectrum across the wavelength axis caused by the redshift. We also see that the 
the lower wavelength portion of the spectra gets more attenuated with increasing redshift. 
This is because the optical depth peaks towards lower wavelengths for this grain size and 
the optical depth increases with redshift. In the second figure we see that with lower dust 
grain size it leads to increasing attenuation. All of this behaviour is explained by the optical 
depth wavelength dependency explained above. 


7.5 Comparing predictions to observed quasar colours 


As illustrated in Figure 10 we are now able to predict the attenuated quasar spectra. To 
compare these model predictions to observations we will be using colours. This is advanta- 
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Figure 17: The redshift distribution of the quasars in SDSS DR7. Here we see that we have 
a larger amount of quasars in the redshift range 0-2, hence we will have a larger statistical 
significance from this redshift range. 


geous since there are a lot of photometric measurements of quasars publicly available (Lyke 
et al. 2020). Additionally since colours are the difference in magnitudes for different filters 
(see Section 4) it is essentially a logarithmic difference between fluxes. This means that any 
systematic effect or factor that is multiplied evenly over the entire spectrum cancels (see 
Equation 14). Also note since our simulated spectrum is in arbitrary flux units by using 
colours such normalization factors will cancel. 


7.5.1 Observed quasar colours 


We will use data from SDSS’s Data Release 7 (DR7) (Schneider et al. 2010), which contains 
105783 confirmed QSO’s with magnitudes in the five SDSS filters (u,g,7r,7, z). Only QSO’s 
with magnitudes over [u, 9, 7,7, 2] = [22.3, 22.6, 22.7, 22.4, 20.5] is considered, corresponding 
to a S/N ratio greater than 5:1 to avoid quasars with too much noise (Richards et al. 2001). 
This cuts the data down to 97153 QSO’s. Quasars which are classified as extended sources 
are also disregarded to minimise the contamination from host galaxies. What is left is 93305 
point sources with sufficient S/N ratio. We also only look in the redshift range 0.5 < z < 3 
since that is the range with the most amount of quasars in the data, resulting in the greatest 
statistical significance. Also below z = 0.5 the data is more likely to contain Seyfert galaxies 
and the quasar spectra at these redshifts may be contaminated by their host galaxies (M6rtsell 
and Goobar 2003; Ostman and Mértsell 2005). The redshift distribution of the SDSS quasars 
is shown in Figure 17. 

To get the colours from the real observations we divide the observations into redshift bins 
of size Az = 0.05 in the redshift range 0.5 < 3. For each bin we calculate, for every filter X 
and Y, the average colour (X — Y)?°° for all of the observed quasars in that bin. Quasars 
that diverge more than 20 from the average colour are rejected to remove contaminations 
from anomalous colours, for example from quasars which has extensive extinction from their 
host galaxy (Mértsell and Goobar 2003). The SDSS flux calibration should take much of 
these effects into account but outliers always exist (Masataka Fukugita 2011). 
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Composite quasar spectrum with SDSS filters 
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Figure 18: The SDSS filter transmission functions superimposed on the composite spectra. 
This shows which part of the spectra is covered by the filter functions. As the spectra is 
redshifted the spectra will move to the right in the plot, changing which sections of the 
spectra are covered by a certain filter. 


However the spectrum blueward of the Lyman-a (Lya) emission line are heavily affected 
by absorption due to the Lya forest, thus decreasing the flux from those wavelengths (Vanden 
Berk et al. 2001). This is apparent when looking at Figure 11 where the error is greatly 
increased blueward of the Lya line. That means at higher redshifts the bluer filters will 
extend into this region making the filters inefficient for those redshifts. To address this we 
only include QSO’s with z < 1.5 for the u-filter and z < 2.3 for the g-filter (Johansson and 
Moértsell 2012). Now we have all 10 different colours for each redshift bin, the result of this 
can be seen in Figure 19. 


7.5.2 Predicted quasar colours 


For the simulated colours we also divide them in the same redshift bins. Then we specify 
the dust grain material and size. For each redshift bin we create a dust attenuated spectrum 
following the steps in Section 7.3. For each redshift bin we then apply the SDSS filter 
functions to the spectrum for that bin to get the simulated fluxes in the SDSS filters. 

We apply the filter functions on the simulated spectrum as: 


a / fas Tr(A)dd (27) 


where f, is the filter flux for the filter x, f, is the composite spectrum flux density and 
Tr(A) is the filter sensitivity depending on wavelength. In Figure 18 we can see the filter 
transmission curves along with the composite spectrum. 

Practically we compute the filter fluxes using the python package Pyphot (Pyphot, M. 
Fouesneau 2016). Pyphot is equipped with tools to calculate and filter flux, and is able to 
account for scenarios where the filter functions extends beyond the available spectrum (as 
can be seen in Figure 18 for the z filter). Once we have the filter fluxes we can easily convert 
these into colours. Again we only calculate colours with the u-filter for bins with z < 1.5, 
and the g-filter for bins with z < 2.3 as explained in the previous section. In Figure 20 we 
see the colours for silicates, the left figure shows all colours for a = 0.1 and the right figure 
shows the i — z colour for several dust sizes. 
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Observed colour over redshift 
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Figure 19: The evolution of the 10 SDSS colours for observations at different redshifts. Here 
we see how the colours change as the quasar spectra gets redshifted, which results in the 
filter functions covering different parts of the spectrum as it gets redshifted (see Figure 18). 
These will later be compared to the expected colours from the model simulations. 
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Figure 20: The colour evolution over redshift. All of the colours for a dust model with a = 0.1 
jm (left) and the (i-z) colour for models with several different dust sizes (right). We see that 
the shape and magnitude of the left figure colours resembles the observed colours of Figure 
19. We also see slight differences in the (i-z) colour between different dust sizes in the right 
figure, especially towards larger redshifts. 


7.5.3 Comparing simulations to observations 


Now that we have all of the colours for both the simulations and the real observed colours 
we can start to compare them. We do this using the following equation 


A(X -Y),=(X -Y)%—-(X-Y)E", (28) 


where (X — Y)°°° is the observed colour and (X — Y)**” is the simulated colour for the filters 
X and Y. This equation will be utilized in the coming section about the statistical analysis 
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Figure 21: Side by side comparison of the observed colours and the colours from the composite 
spectrum without the dust model applied to the spectrum. Here we see the colours are fairly 
similar both in terms of shape and magnitude, which is expected since they are both the 
quasar colours. 
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Figure 22: The colour evolution over redshift for the (i-z) colour. The left figure shows the 
observed colours (dashed) and the colours of the composite spectrum without the dust model 
applied to it (full). The right figure shows the colour for a couple of different dust models 
being applied to the composite spectrum, as well as the observed colour. 


of different dust models. In Figures 21, 22 and 23 a couple of plots can be seen showing 
different colours and comparing simulated colours to observed colours. 

As can be seen in Figures 22 and 23 below about redshift z = 1.3 we see that the observed 
colours are bluer compared to the simulated colours. This is an error in the analysis. Due 
to this offset any dust attenuation applied to the simulated colours below this redshift will 
lead to a negative correction since dust attenuation leads to reddening. This is likely an 
error caused by the composite spectrum. As explained in Section 7.2.2 there are possible 
contaminations from intergalactic dust, as well as other sources of errors. These effects could 
lead to this offset, possibly producing a bias in the result. As will later be seen in Section 
7.8 the error from the composite spectrum is the leading error in the analysis, i.e reducing 
the errors in the composite spectrum is of great importance to improve on the result. 
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Figure 23: The colour evolution over redshift for the (r-i) and (r-z) colours. The figures shows 
the colour for a couple of different dust models being applied to the composite spectrum, as 
well as the observed colour. 


7.6 Testing colour calculations with varying degree of complexity 


We tested the implementation of the quasar colour calculation by using different spectra 
which had varying degrees of complexity in the dust model which were applied to them. 
We then calculated the colours of the spectra and analysed the result. These spectra are as 
following: 


1. A baseline spectrum with constant flux, basically a horizontal line. 
2. The composite spectrum, without any dust attenuation at all. 


3. The composite spectrum with dust attenuation but with constant interaction efficiency, 
leading to a constant cross section for a certain dust grain size. 


4. The composite spectrum with changing sigma, the model we are looking into. 


We chose these different test scenarios to highlight different aspects of the physics model. 
The aim of this section is essentially to investigate the colour analysis and see if it agrees 
with theories and equations. This is useful in the search for errors in the calculations, which 
would be found if the actual results would differ from the expectations from theory. For 
example we want to test that the shape of the colour curve comes from the change in flux 
over the spectrum, that the wavelength dependency in the optical depth comes from the 
extinction efficiencies etc. We do this by starting from a very simplified model, the baseline 
spectrum, which is not a real physical model but just constructed for testing purposes. Then 
add complexity, like using the composite spectrum or adding the dust attenuation model. 
For each successive step we analyse the results to see whether the produced results agree or 
not with the expectations, and if the results make sense in the context of what we are doing. 
For these tests the model which assumes a constant amount of dust grains was used (e.g 
Model 1, see Section 7.3.1) since this model produces the most significant changes in colours 
depending on the dust size. Note that this is the only section using this model. All of the 
plots are for silicates and at redshift z = 2. Below we further explain each spectrum as well 
as the results from the colour calculations. 


7.6.1 1. Baseline spectrum 


We assume the spectrum to just be constant over wavelength, which also means the spectrum 
is invariant under redshift. Thus the colour we calculate will stay the same indifferent of 
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redshift resulting in constant horizontal lines in the colour over redshift plot, which is exactly 
what we see in Figure 24. This is not a physical model, just used to test the technical 
implementation of the colour calculations. 


7.6.2 2. Composite spectrum with no dust attenuation 


To understand the impact of the spectral features, we use the composite spectrum without 
any dust model applied to it, as described in Section 7.2.1. This means the that the colour 
calculations will vary over redshift as the spectrum is no longer constant over wavelength, 
meaning as we shift the spectrum in the wavelength dimension the colours will change as the 
flux changes. Thus we expect the colour calculation of the quasar composite spectrum to 
closely resemble the observed quasar colours from SDSS, which is what we see in Figures 24 
and 22(left) where both the magnitude and shape of the colour curve resembles the observed. 
This is a good confirmation that at least part of the colour calculation and its application 
on to the composite spectra are correct. 


7.6.3. 3. Composite spectrum with dust attenuation with constant cross-section 


As a next step we introduce our dust attenuation model, but for simplicity we first set the 
interaction efficiency to a constant. This results in a constant cross-section since we calculate 
for a certain dust grain size (a = 0.1 um). Thus we remove the wavelength dependence from 
the optical depth (equation 22) and our model which results in the dust attenuation just 
being a constant which is applied equally to the entire spectrum. When we calculate the 
colours, which is a logarithmic difference between parts of the spectrum, this constant will 
get removed (see Equation 14). The result from this should be that the colours are identical 
to the case where we have no dust attenuation, which is exactly what is seen in Figure 24, 
where the colours perfectly overlap. If they would not overlap it would mean we had some 
wavelength dependence in our optical depth outside of the cross section, which we should 
not have. Yet another promising result. 


7.6.4 4. Composite spectrum with dust attenuation 


Finally we use the full dust model which has a wavelength dependence. Here we are calcu- 
lating the colours for two different dust grain sizes. Generally for larger optical depth we 
are expecting the colours to be further reddened by the intergalactic dust, e.g the colours 
are larger for a larger optical depth. A general feature of the model we are looking at in 
this section (Model 1) is that with increasing dust grain size the optical depth will increase 
(see Figure 27). This is because we assume a constant amount of dust grains, thus with 
increasing dust grain size we will effectively get more dust in the universe (in mass units). 
The result of this is that larger grain sizes will lead to larger colours. Another thing affecting 
the result of the colour calculation is how the optical depth changes over redshift. Since the 
light from the object has to travel through more dust the further away it is the optical depth 
increases with redshift (again see Figure 27). These features leads to us expecting the dust 
attenuated colours to be slightly (or not so slightly depending on dust size) larger than the 
colours from the composite spectrum, with the reddening increasing towards larger redshifts. 
This is exactly what we see in Figure 24 where the spectra with our dust model are more 
reddened, with larger reddening for larger dust grain sizes. Again a good confirmation that 
the colour calculations work and does what we expect when we apply our dust model. 
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Figure 24: Comparing models with different complexity against observations using the 
colours r-z (left) and i-z (right). We see the baseline as a horizontal line, the spectra without 
dust and the one with constant o overlaps, and the spectra with grain sizes has increased 
reddening with larger grain sizes. We also see the spectra with a = 0.05 zm are comparable 
to the observed colours. 


7.6.5 Conclusion from the testing 


The results from the tests of the different models with varying complexity confirms that the 
model correctly accounts for the relevant effects. The results shows: 


e The colour calculation depends on the changing flux over wavelength of the spectrum, 
as expected. 


e The colours of the composite spectrum matches the observed SDSS colours in both 
shape and magnitude. 


e Our only wavelength dependency in the dust model comes from the extinction efficiency 
from the cross section, as it should. 


e Applying the dust model from Model 1 to the composite spectrum reddens the colours, 
with further reddening for larger dust grain sizes and towards larger redshifts. 


7.7 Discussion of colour comparison 


Now lets look at the results of the colour analysis. The shape of the colour vs redshift plots 
are caused by the shift of the spectrum along the wavelength axis as the quasar is redshifted. 
Since the wavelength span of the filter functions is set, as the quasar spectra are redshifted 
different part of the spectra will move in and out of the filter functions. This then leads to the 
magnitudes, and consequently the colours, changing over redshift. For example the spectra 
might peak in one particular filter for a particular redshift. But as the spectra is shifted 
towards different wavelengths the peak may move into the wavelength span of different filter 
functions. 

In Figure 20 we see all of the simulated colours and the simulated (i-z) colour for a 
couple of different dust sizes. In the right figure we see that the colour varies slightly for 
each different dust grain size. This is because of the different wavelength dependency in the 
optical depth for different dust grain sizes as seen in the right figure in Figure 14. The dust 
grain size that leads to the most dust attenuation, i.e largest optical depth, over a certain 
filter will cause the most reddening and thus have the largest colour. 
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We can now look at the comparison between the observed and simulated colours. In 
Figure 21 we see a side by side comparison of the observed colours vs the simulated colours 
for only the composite spectrum, without our dust model applied to it. We see that the shape 
of the colour curve matches fairly well. Then in Figures 22 and 23 we see a comparison of the 
(i-z), (r-i) and (r-z) colours. In Figure 22 we firstly see the simulated composite colours and 
the observed colours compared to each other in the left figure, and in the right we compare 
the observed and simulated colour for different dust grain size models. 

We see that certain dust size models produce colours that better match the observed 
colours than others. This is what we look at to determine which dust size model is the best fit 
to the real observations. As previously discussed we also see that all of the simulated colours 
are always more reddened towards lower redshifts below about z = 1.2. This points toward 
some possible systematic error in the composite spectrum towards lower redshifts. It might 
be an offset created by the effects of intergalactic dust on the generation of the composite (see 
Section 7.2.2). This offset might cause a bias in the results and to produce more definitive 
results improvements are needed taking into account or removing the systematic effects in 
the composite spectrum. 


7.8 Analyzing parameters to fit simulations to observations 


We are now ready to look at different dust models, i.e different dust grain sizes and compo- 
sition, and see how they compare to observations. We do this by performing a Pearson’s y? 
test (Pearson 1900). Essentially what we do is calculating the square difference between all 
of the observed and simulated colours for each redshift bin for each dust model. This is done 
by the following y? equation (Johansson and Mértsell 2012): 


N 
P= >> AK -Y)VG' A(X -Y);, (29) 


ij=l 


where N is the total number of redshift bins, and 7 and j are the different bins. V;,; is the 
covariance matrix and 
2 2 2 
Vig = Feali + F comp + bij obs5 (30) 


where O¢qi;j is the calibration uncertainty of the photometric system and is given as ~ 0.02 
(M. Fukugita et al. 1996). comp is the uncertainty of the power-law fits to the composite 
spectrum and is given as ~ 0.10 in Vanden Berk et al. 2001, and lastly oo», is the uncertainty 
of the SDSS DR7 colours and is given to be ~ 0.003 in Schneider et al. 2010. 

The idea for this statistical analysis is that the dust model which best fits to the obser- 
vations would produce the smallest +7. Hence we can look at x? to determine how good the 
different models fit the observations and thus constrain what are the more likely dust sizes 
for each dust composition. In Figure 25 we see the result from this analysis for both silicate 
and graphite for dust models with Qaus¢ = 5 x 107°. For these dust models the results favors 
dust grains around size a = 0.05 um for graphite and a = 0.08 ym for silicate. Also silicate 
produces a smaller x? meaning that silicate produces the better fit to observations. 

We also performed a 2d grid search for both the grain size and Qgus: for each dust model 
for silicate. The result of this can be seen in Figure 26 which shows the logarithmic value of 
the y? as the colour with smaller values being darker, i.e the darker the colour the better the 
fit to observations. This calculation show a favor for dust grain sizes around a = 0.08 zm and 
Qaust around 5— 7 x 1076. This result is in agreement with previous estimates (Johansson 
and Mortsell 2012). 
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Figure 25: The results from the x? calculations for both graphite and silicate. This is a 
measure of how closely the dust model matches observations where the smaller x? the better 
that model fits to observations. This means that the grain size that produces the minimum 
in y’ is the dust size that produces the colours that is the closest to the observed data. 
The reason behind the "step-like" curve is the extinction efficiency data. This data is only 
available for certain dust sizes leading to the "step-like" shape. 
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Figure 26: The two dimensional colour plot calculation spanning over both Qaus; and the 
dust grain size for silicate. The colour shows the logarithmic y? value for that grid, i.e the 
darker the colour the more probable the dust model is and the minimum is marked with a 
white dot. The colour bar shows the logarithmic y? value. 
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8 Discussion of the statistical analysis and findings 


Now we come to the statistical analysis part of this thesis. Here we compute, for each dust 
mode, the total difference between the simulated and observed colours for all of the colours 
for all of the bins and sum them up. This was done for several dust models with varying dust 
grain size, composition and Qgys¢. The model that produces the smallest difference is then 
the model that is best fitted to the observations and the dust model that best represents 
real intergalactic dust. In this process we take into account the errors from the composite 
spectrum, the colour calculations and the error in the observed data. As seen in Section 7.8 
the combined error of this calculation is determined by the error of the composite spectrum, 
and the error from the observed data and the colour calculations are comparably insignificant. 
This means we must first reduce the errors in the composite spectrum before reducing the 
other errors would contribute to significant reductions in the total error. Thus using a newer 
data set with improved errors would not contribute much in the improvement of the results. 
To produce more definitive results improvements of the composite spectrum are required, 
specifically taking into account or removing systematic effects such as the possible effect of 
intergalactic dust in the composite spectrum. 

The result of the statistical analysis suggests that a dust model with Qyys¢ = 5 x 107° 
favors dust grains around size a = 0.05 wm for graphite and a = 0.08 ym for silicate. Both of 
these results are within the limits given in Johansson and Mortsell 2012, with the graphite 
size being just above the lower limit. This analysis also suggests that dust models with 
silicates are favored, as the y? was smaller for silicate compare to graphite. When the 
statistical analysis was done over both ranging dust sizes and OQdust the results suggests that 
dust models with the same dust sizes as mentioned above with a Qq,s; around 5 — 7 x 107° 
are favored. This is also in rough agreement with previous estimates of Quus: (Masataka 
Fukugita 2011; Johansson and Mortsell 2012). This shows that intergalactic dust exist, and 
has a notable effect on cosmological measurements and favors dust grains consisting of silicate 
with a size around 0.1 um. 


9 Conclusion & summary 


In the current cosmological paradigm we have a universe which is expanding at an accelerated 
rate. This expansion is supported by the observations of Type Ia supernovae. However in 
order to get accurate measurements one has to take into account the systematic effects 
and errors that may influence the observations. The reddening and dimming caused by 
the interaction between the light emitted by the supernovae and the intergalactic dust on 
the way is one of those systematic effects. As the precision of cosmological measurements 
increase the effect of such systematic effects becomes more important to take into account. 
Currently there are also a tension in the value of the cosmological parameter Hp, as two 
different methods produce different results. Systematic effects caused by intergalactic dust 
could be contributing to this tension (Mortsell et al. 2021) and understanding these effects 
could help in understanding this tension. 

In this thesis we analysed the influence of intergalactic dust on cosmological measure- 
ments. To do this we constructed a forward model depicted in Figure 10 where we had a 
simulated quasar emit light through the intergalactic medium where it encountered inter- 
galactic dust. For the quasar spectra we used a composite spectrum (see Section 7.2.1) and 
how the light interacted with the intergalactic dust grains is governed by the equation shown 
in Section 5.5. We were thus able to apply these equations to the simulated quasar spec- 
tra for a chosen redshift and dust model to create the intergalactic dust attenuated quasar 
spectra (see Section 7.3). For the different dust models we assumed spherical dust grains of 
either silicate or graphite, and could change the dust size, dust composition and OQ, for 
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each simulated spectra (see Sections 5.5.1 and 5.5.2). Once these algorithms were in place we 
could compare these simulated spectra to real observational data. This was done using SDSS 
observations and comparing the simulations to the observed data using colours (see Section 
7.5). Lastly we did a statistical analysis for different dust models, essentially measuring the 
difference between the simulated colours and the colours from the observational data (see 
Section 7.8). The dust model that produced the smallest difference was thus the dust model 
that best fitted the observed data and were the dust model best portraying real intergalactic 
dust. To produce more definitive results improvements of the composite spectrum used in 
the analysis are required. Specifically surrounding the possible systematic effects. 

The result of this analysis favors grains consisting of silicate with a dust grain radius 
around 0.1 ym and Qgy.; around 5 — 7 x 107°, which is in agreement with previous results 
(Johansson and Mortsell 2012). We can therefore conclude that intergalactic dust exist and 
they have an effect on measurements. 


10 Outlook & prospects 


This work did some assumptions, and improvements on the method is possible. Since this 
thesis assumed equal number density of dust in all directions one could introduce a inho- 
mogeneous distribution of dust grains mapping the distribution of galaxies to get a more 
accurate dust distribution model. One could use this to construct composite spectra with 
less possible effects from intergalactic dust, as well as taking quasars from lower redshift to 
construct the composites (Ostman and Mértsell 2005). It is also possible to use this distri- 
bution of dust to produce better intergalactic dust models and more accurate corrections for 
the effects of intergalactic dust on observations. One could also use a varying composition of 
different materials and grain sizes, where this thesis only assumed a certain grain material 
and dust grain size. There is also the possibility of further analysing the spectral features 
intergalactic dust causes, and look at something else besides colour to compare simulations 
to observations. Lastly one could also investigate the errors further, how large are the sys- 
tematic effects in the composite spectrum or the SDSS data. As stated in previous sections 
the largest improvements of the results would come from improvements of the composite 
spectrum used. It would thus be fruitful to do this analysis with an improved composite 
spectrum. 
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A Results of Model 1 


All of the plots shown below are the results from Section 7 but with Model 1 instead of 
Model 2. 
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Figure 27: The optical depth as a function of redshift, r(z). First plot with constant wave- 
length and varying dust size and the second plot with constant dust size and varying wave- 
length. Both the dust grain size and wavelength are in wm units. 
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Figure 28: How the optical depth depends on wavelength with varying redshift for silicates. 
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Figure 29: How the optical depth depends on wavelength with varying dust size in wm for 
silicates. 
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Composite spectrum with dust attenuation (z=0.5) Composite spectrum with dust attenuation (z=1) 
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Figure 30: Composite spectra with dust attenuation caused by silicates for varying redshifts. 
All of the spectra has a dust size of dgust = 0-luum. 
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Composite spectrum with dust attenuation (a=0.01 um) 
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Figure 31: Composite spectra with dust attenuation caused by silicates for varying dust sizes. 


All of the spectra has a redshift z = 2. 
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Composite spectrum with dust attenuation (a=0.05 um) 
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Figure 32: The results from the x? calculations for both graphite and silicate. This is a 


measure of the how closely the dust model matches observations where the smaller \? the 
better that model fits to observations. 
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Observed colour over redshift 
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trum without dust attenuation. 
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Figure 33: Several plots showing the result of the colour calculations, along with some 
comparisons between the observed and the simulated for the composite spectrum with no 


dust attenuation.) 
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